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We show by using the method of matched asymptotic expansions that a sufficient condition can 
be derived which determines when a local experiment will detect the cosmological variation of a 
scalar field which is driving the spacetime variation of a supposed constant of Nature. We extend 
our earlier analyses of this problem by including the possibility that the local region is undergoing 
collapse inside a virialised structure, like a galaxy or galaxy cluster. We show by direct calculation 
that the sufficient condition is met to high precision in our own local region and we can therefore 
legitimately use local observations to place constraints upon the variation of "constants" of Nature 
on cosmological scales. 
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OA , I. INTRODUCTION 

1 . The last few years have seen a resurgence of widespread interest in the possibility that some or all of the fundamental 

'constants' of Nature might be varying over cosmological timescales. To a large extent the revival in this field has come 
, about because of the need to explain and understand the studies of relativistic fine structure in the absorption lines 
formed in dust clouds around quasars carried out by Webb et al., Using a new "many-multiplet" method which 
, exploits the information in many wavelength separations of absorption lines with different relativistic contributions 
J> ■ to their fine structure considerable gains in statistical significance were achieved. From a data set of 128 objects at 
] redshifts 0.5 < z < 3, Webb et al. found their absorption spectra were consistent with a relative in the value of the 
, fine structure constant, aem{z), between those redshifts and the present of Aaem/ctem = ctemiz) — aem(0)/aem(0) = 
— 0.57± 0.10 X 10~^. In the seven years since their results were first announced, extensive analysis has yet to identity 
any systematic effect that could explain either its magnitude or sign. However, a small study of 23 absorption 
\^ ' systems in one hemisphere between 0.4 < z < 2.3 by Chand et al., found a result consistent with no variation: 
[ Aaem/aem — —0.6 ± 0.6 X 10^^ but this study uses only a simplified version of the many-multiplet method in its 
' analysis and concerns remain about calibrations and the noisiness of the data fits. Later this year, a major effort to 
^H, produce a very large new data set should be reported and will clarify the status of these earlier investigations. In the 
I ■ meantime there have been other astronomical checks on the constancy of aem, using a variety of different techniques; 
^jpj' however, they have yet to reach the accuracy achieved using the many-multiplet method. All found results 's'l that 
• • : are consistent with no variation but none could have seen the variation reported by Webb et al. If aem can, and does, 
' vary in time then it would seem natural to suspect that some other fundamental 'constants' do so also. Recently, 
k> using a study of the vibrational levels of H2 in the absorption spectra of quasars, Reinhold et al., 0, reported a 3.5a 
5_] indiction of a variation in the electron-proton mass ratio, /i = nie/mp over the last 12 Gyrs: A/x//x = 2.0 ±0.6 x 10~^. 
[ This result combines high-quality astronomical data with improvements in the measurement of crucial laboratory 
wavelengths to new levels of precision. It reports a variation at a level comparable to that claimed for aem (which 
partly reflects state of the art precision in spectroscopic measurements) but is theoretically much harder to understand. 
Variations in aem and /i both lead to violations of the Weak Equivalence Principle (WEP) and these are expected 
to be unacceptably large if Afi/fj, varies at the 10~^ level. Typically, the experimental constraints on WEP violation 
from direct experimental studies of freefall and the study of the relative motions of the Earth- Moon system lead to 
upper bounds on the relative differential accelerations of different materials which are O(10~^^). As explained in the 
detailed study by Barrow and Magueijo 0, the Webb et al observations of varying aem predict violations at the 10"^"^ 
level but the Reinhold et al. observations of varying fj, lead us to expect violations at the 10~^ level. 

In any study of varying constants, the data used to constrain our theories which allow variations to occur comes 
from a number of very different environments and scales: with densities differing by a factor of 10"^*^ or more, and 
spanning some 12 billion years. In order to be able to use all of the information available we need to know how 
the results of local laboratory experiments, terrestrial or solar-system bounds from the Oklo natural reactor, and 
from isotope ratios in meteorites, are related to data coming from astronomical observations on extragalactic scales. 
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This is the 'Local vs. Global' problem for varying-constants. It is an important problem and yet most commentators 
invariably assume that the local and cosmological observations are directly comparable [IQiJ . This is strong assumption 
and is almost invariably made without any proof. A priori, it is not obvious that this assumption is true; indeed, in 
many other theories, not least that for gravity itself, it is not: we do not expect to be able to measure the expansion 
of the Universe by observing an expansion of the Earth. In this paper we describe the first rigorous proof of why, 
in almost all varying-constant theories, local experiments will also 'see' any variations in 'constants' which occur on 
cosmological scales. 

II. GENERAL THEORY 

Before we can solve the 'Local vs. Global' problem, we need to introduce the general way in which a constant is 
promoted to become a dynamical quantity consistent with Einstein's conception of gravity. In general, a constant, 
C, is allowed to vary by associating it with some scalar field or "dilaton", (p, i.e. C — > C{(j)). We usually assume 
that the scalar field theory associated with has a canonical kinetic structure and the variations of this scalar field 
contribute to the spacetime curvature like all other forms of mass-energy. Variations in (j) must also conserve energy 
and momentum and so their dynamics are constrained by a non-linear wave equation of the form 

where is associated with the variation of one or more 'constants', Cj, via a relation Cj = fj.ipi't') — '^fj{4')/d4'- 

The Lj(ek,Pk) are some linear combinations of the density, Sk, and pressure, pk, of the fc*'* species of matter that 
couples to the field 0. Included in this formulation are all standard theories for varying constants, like those for the 
variation of the Newtonian gravitation 'constant' G, aem, and the electron-proton mass ratio, as described in refs. 

mm 

For local observations to be directly comparable with cosmological ones we need to know the conditions under 
which 

(j){x, t) K (j)^{t) 

to some specified precision, with x taking values in the solar system, where the subscript c labels the large-scale 
cosmological value of the field (j). The validity of this approximate equality, the accuracy to which it holds, and the 
accompanying preconditions needed to support its validity are the subject of the rest of this paper. 

Prior to the onset of the matter era the universe is homogeneous to a very high precision inside the horizon. Any 
inhomogeneities that do exist, and the evolution of within them, can be consistently and accurately described by 
linear perturbation theory and (/) ~ 0c holds. But the study of the evolution of "constants" becomes mathematically 
challenging when linear theory breaks down and the inhomogeneities become non-linear. This only starts to occur 
during the matter era; at these epochs it is an acceptable approximation to consider the Universe to be comprised 
of only pressureless dust (baryonic and dark matter), density e, and some cosmological constant, A. We usually 
expect that the scalar field will couple only to some fraction of the total dust density; for example, in varying-agm 
theories it couples to the fraction that feels the electromagnetic force, and in varying theories it couples only to 
the electron density We will assume, as is almost always the case, that the fraction of matter to which it couples 
is approximately constant during the epoch of interest. Under these simplifications eqn. ^reduces to: 

where k = Stt, c — G = 1, V{4i) is the dilaton potential, and B^^{(j)) is the effective dilaton-to- matter coupling. In 
what follows we shall assume that the varying-constant evolves according to the above conservation equation, which 
is certainly true of Brans-Dicke theory, BSBM and BM varying-/! theory. In what follows we will further assume that 
the cosmological value of 0, denoted by (f>c, is sufficiently far away from any extrema of the matter coupling, B{(j)). 
We also demand that V^{(p) not be too large. Our conditions on B{(p) and V{(p) are summarised as follows: 

for all values of (f> within the range of interest (i.e. those that can be reached from the evolution of some given 
initial data). The condition on the matter coupling is usually equivalent to \B^^^{(j)c)\ <C 1 and B^^{(j)c) 7^ 0. The 
condition on the potential must hold for = 0c to prevent the varying "constant" evolving at an unphysically fast 
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rate cosmologically; the assumption that holds everywhere will then be valid provided that V4,{4>) is suitably flat. As a 
result of this final assumption our results will not apply to Chameleon field theories, ^ 12j . With the major exception of 
theories with "Chameleonic" behavior, our model includes almost all physically viable proposals for varying-constant 
theories. Our results are also applicable to any scalar-field theory, not just those that describe varying-constants, 
provided that the scalar satisfies a conservation equation of the above form. For a general matter distribution the 
dilaton conservation equation is a second-order, non-linear PDE, and there is no reason to suspect that it should be 
easily solvable, or indeed analytically solvable at all. Even numerical calculations will generally be difficult to set-up 
and control. Cosmologically, we assume homogeneity and isotropy which leads us to a FRW background and to a 
solution for — (pdt). Under these specifications the conservation equation for (j) reduces to an ODE in time, and 
can be solved. The other scenario in which it reduces to an ODE is near a spherically symmetric, static body which 
couples to the dilaton strongly enough so that any temporal gradients of (f) are negligible compared to the spatial 
ones. In these cases we can easily find the leading-order static mode of cj), but to find the temporal evolution of (f> 
we need to enforce the boundary condition that it match up to its cosmological value at large distances. The central 
technical problem is that the local, static solution was found under the assumption that the temporal derivatives were 
negligible, and at infinity this is no longer the case. Indeed, to get to a region of space where we know (j) « (t>c{t) we 
must certainly pass through some zone where the temporal and spatial gradients of are of comparable magnitude. 
As soon as we reach this zone, the assumptions under which the local solution was derived break down. In short: we 
cannot consistently apply the boundary condition at infinity to the approximate local solution since spatial infinity 
is far outside the range of validity of that approximation. We shall express this idea more formally below and see 
that it is associated with the fact that the local asymptotic approximation is not uniformly valid. To circumvent this 
problem, created by the presence of multiple length scales, we will use the method of matched asymptotic expansions. 



III. MATCHED ASYMPTOTIC EXPANSIONS 



Second-order, non-linear PDEs are difficult (and often impossible) to solve exactly. However, if one can identity 
some small parameter, S, in the problem then it is usually possible to find an expansion in d which is formally 
asymptotic to the solution in the vicinity of some fixed point. An approximation /„(a;)7„((5) is asymptotic to a 
fimction f{x,6) as (5 ^ iff 

fM{x)jMiS) 

for fixed x. If this definition holds for all M then we write: 



as J ^ 0, 



f{x,S) - ^/„(a;)7„(a;), 

and fn{x)^n{x) is an asymptotic expansion of f{x, (5) as 5 — > for fixed x. The sum here is a formal sum, since in 
general it will not converge; however, as a result of the defining property of asymptotic expansions we will in general 
need only the first few terms of the sum to obtain a very good approximation to /(x, S) at x. Asymptotic expansions 
are unique for each x, but it is also clear that an expansion that is asymptotic to f{x, S) for some range of x, with 
X ^ 0{l) say, will not in general be valid in some other range of x, usually x ~ 0(1/(5) or a; ~ C>{S). In these cases the 
expansion is said to be not uniformly valid. If an expansion has arisen as an approximation to the solution of a PDE 
and is not uniformly valid, then the PDE is said to exhibit singular behaviour. Such behaviour is often associated 
with the presence of two or more very different length or times scales in the problem. This is precisely the case in the 
'Local vs. Global' cosmological problem, where the length scale of the local inhomogeneity is very much smaller than 
the Hubble scale which defines the cosmological background. 

We can proceed with such problems by constructing two (or more) asymptotic approximations to the solutions 
which are valid for different ranges of x, e.g. for x ~ 0{1) and x/S = (, ^ C(l); with 

f{x,S) r^Y.n=ofn{x)^n as (5 0, x fixcd, (2) 
f{x, S) - Y.n=o 9n{i)5n as (5 ^ 0, ^ = x/5 fixed, (3) 

and solving the PDE order by order in 5 for both expansions w.r.t. to some boundary conditions. We call expansion 
P|) the outer solution, and ^ the inner solution. The inner expansion is not uniformly valid in the region ^ = Oil/ 5), 
as the outer one is not valid where x = 0{5). Because of these restrictions on the size of x, we will only be able 
to apply a subset of the boundary conditions to each expansion; in general, we will therefore be left with unknown 
coefficients in our asymptotic approximations. This ambiguity can be lifted if there is some intermediate region. 
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e.g. X ~ 0{S^^^) where they are both vaHd, by appealing to the uniqueness of asymptotic expansions and matching 
the inner and outer solutions there. In this way we can effectively apply all boundary conditions to both solutions. 
This is the method of matched asymptotic expansions (MAEs) . Its application to problems in general relativity was 
pioneered by Burke, Thorne and D'Eath in the 1970s. For a fuller account of MAEs we refer the reader to refs. 
and [Hi. 

IV. GEOMETRICAL SET-UP 

The experimental bounds on the permitted level of violations of the WEP due to the presence of light scalar fields 
demand that the dilaton field couples to matter much less strongly than gravity, so 

|S 0I « 1. 

As a result, the dilaton field is only weakly coupled to gravity, and so its energy density and motion create metric 
perturbations which have a negligible effect on the expansion of the background universe. This feature allows us to 
consider the dilaton evolution on a fixed background spacetime. In this work we go further than we did in ref. [Tll | 
and consider not only the extent to which condition s:! 0c is satisfied near the surface of some spherical virialised 
over-density of matter, e.g. a the Earth, a star, black-hole, galaxy or galaxy cluster, but also the degree to which it 
is valid during the collapse of an over-dense region. We will, however, treat the two cases separately. 

In the first case, we shall refer to the virialised over-density as our 'star' and take it to have mass m and radius Rs 
at some time of interest t = Iq. Although we require that the 'star' itself be spherical, we do not demand that the 
background spacetime possess any symmetries. We do require however that, at i = to, the metric is approximately 
Schwarzschild, with mass m, inside some closed region of spacetime bounded by a surface at r = i?s; this region is 
called the interior. The metric for r < Rg is left unspecified. We allow for the possibility that r — Rg is a. black-hole 
horizon. 

In the second case, we only consider the case where the spacetime is spherically-symmetric, label the mass of the 
collapsing region by m, and assume that its spatial extent is small compared to the Hubble scale. We also demand 
there are no black-hole horizons in the interior of the collapsing region. By using the results of the first case, however, 
we can, in some cases, allow for the formation of a horizon. 

In both cases we demand that: 

• Asymptotically, the metric must approach FRW and the whole spacetime should tend to the FRW metric in 
the limit m —> 0. 

• The spacetime is approximately FRW in some open region that extends to spatial infinity, this is called the 
exterior. 

We are concerned with spacetimes where the matter is a pressureless dust of density £, with cosmological constant, 
A. We further require that the motion of the dust particles be geodesic. In the spherically-symmetric case, all such 
solutions to Einstein's equations with matter fall into the Tolman-Bondi class of metrics (for a review of these and 
other inhomogeneous spherically symmetric metrics see ref 18]), however when condition of spherical symmetry is 
dropped, the general solution is not known. We can simplify our analysis greatly, however, we specify four further 
requirements: 

1. The flow-lines of the background matter are non- rotating. This implies that the flow-lines are orthogonal to a 
family of spacelike hypersurfaces, St. 

2. Each of the surfaces St is conformally flat. 

3. The Ricci tensor for the hypersurfaces St, ^'^^i?a6, has two equal eigenvalues. 

4. The shear tensor, as defined for the pressureless dust background, has two equal eigenvalues. 

These conditions are automatic if spherical symmetry is required, and in general they specify the Szekeres-Szafron 
class of solutions, 14, 15j, of which the Tolman-Bondi solutions, 20, 2lj, are the spherically symmetric limit. We 
require the inhomogeneity to be of finite spatial extent, this limits us to consider only the quasi-spherical Szekeres 
solutions, which are described by the metric: 

„ „ (l + iy rR)^ i?^„dr^ o o / o ON 

ds2 = dt^ - ; , ''^ R^e^'' (dx^ + dy^) , 4 

1 — k(r) ^ ' 
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where := v_r/R.r and 

e""^ = A{r){x^ + y^) + 2Bi{r)x + 2S2(r)2/ + C(r), 



and: 



i?2 = -fc(r) + 2M{r)/R + \kR^ 



In this quasi-spherically symmetric subcase of the Szekeres-Szafron spacetimes the surfaces of constant curvature, 
(i, r) = const, are 2-spheres |22; however, they are not necessarily concentric. These 2-spheres have surface area 
AttR'^, and so we deem 7? to be the physical radial coordinate. In the limit u^r — *■ 0, the (i, r) = const spheres becomes 
concentric. We can make one further coordinate transformation so that the metric on the surfaces of constant 
curvature, {t,r) = const, is the canonical metric on i.e. d^^ +sin^6'd0^: 

X = 2{A{r)x + Bi{r)), 
y^Y ^ 2{A{r)y + B2{r)), 



where X + iY ^ e"^ cot 6/2. This yields 

X2 + y2 + 1 



X^{X^ +Y^ -l) + 2X^X + 2\yY 
r\x.y = v2 i I 1 ^ ^z[r) COS 9 + Xx[r) siii9costp + Xy{r) sin6'sin(/j. 



where we have defined: 



A.(r) :- -J, A.(r) := (^-^ ) ^' ^^(^) ^= j 

With this choice of coordinates, the local energy density of the dust separates uniquely into a spherical symmetric 
part. Eg, and and a non-spherical part. Ens'- 

£ = es{t, R) + £ns{t, R, 0, Lp), 

where: 



1 + V 



,R 



We define M = M r/R.r- We use the remaining freedom to choose r to demand that r ^ R sd t = t^. 

In the virialised case, we follow the conventions of our earlier papers and write M{r) :— m + Z{r). where m is the 
gravitational mass of our 'star'. In the spherically-symmetric = and the metric is of Tolman-Bondi form: 

1 — k{r) ^ ' 

V. VIRIALISED CASE 

In ref. pTl | we considered whether <j){x,t) « 4>c{t) near the surface of some virialised over-density of matter, which 
might be a planet, a black-hole, star, or cluster of galaxies. In general, to specify initial data for the Szekeres-Szafron 
solution we must give both the energy density on some initial hypersurface, ke, and the spatial curvature of that 
hypersurface (given by k{r)). In we considered the two sub-cases of the full Szekeres-Szafron metric, compatible 
with our geometric set-up, where the solution is completely specified by giving the energy density, k£. : 
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• The 'Gautreau' case: the hypersurfaces t = const are spatially flat, fc(r) = 0. In this case the big-bang is not 
simultaneous along the past world lines of all geodesic observers; i.e. it does not occur everywhere for a single 
value of i, at i = say. In these cases the flow lines of matter move out of our 'star' - and the mass of 'star' 
decreases of time. 

• The simultaneous big-bang case: the big-bang singularity occurs at t = for all geodesic observers. In these 
cases the flow lines of matter move into the 'star' - and its mass increases with time. 

The first of these cases is simpler to analyse but the second is more physically reasonable, since we expect the gravity 
to pull matter onto our star rather than expel it. For this reason we will only explicitly consider the simultaneous 
big-bang case in this paper. The results for the Gautreau case are very similar and the simultaneity of the big bang 
is not a significant factor for the late-time evolutionary problem that we are considering. 

The 'interior region', which is immediately outside the surface of the 'star', is approximately Schwarzschild, and so 
an intrinsic interior length scale, Lj, of a sphere centred on the Schwarzschild mass with surface area 47ri?g is defined 
by the Riemann invariant: 

1/4 R^^^ 
' {2mf'^ 

In the asymptotically FRW, or exterior, region, the intrinsic length scale is proportional to the inverse root of the 
local energy density: I/^/kEc + A, where Ec is the total cosmological energy density of matter. We shall assume that 
the FRW region is approximately flat {k = 0), and we define a length scale appropriate for this exterior region at 
epoch at t = to equal to the Hubble radius, which is defined by the inverse Hubble parameter at that time: 

Le = 1/Ho. 

For realistic models Le ^ Lj and so we define (5 to be a small parameter given by 

6 = Li/Le. 

We assume that the whole spacetime metric is Szekeres-Szafron and define dimensionless coordinates appropriate 
to both the interior and exterior near some epoch of interest at time t = t^. In the interior: 

T = L-j\t-ta), ^ = R-'R, 

where R is the physical radial coordinate; T and ^ are 0(1) in the interior and we take the ratio 2m/ Rg to be fixed. 
It is also helpful to define 

rj = (^3/2 _ 3T/2y^^ ; RsV = r + 0(5", S"^/^). 

In the exterior we define: 

T = Hot, p = HQr, 

where r is the unphysical radial labelling coordinate used in the metric 0] We define the interior limit by (5 — > for 
fixed T and ^, and the exterior limit by i5 ^ with r and p held fixed. 

A. The Exterior Limit 

According to our prescription that the metric be FRW to zeroth order in 6, we write 

and 



H^^\,^6n,{p)+oiSn: 
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where s and p are positive numbers which depend on the particular form of the initial matter distribution. The 
exterior expansion of k(r) can be found using the exact solutions for the Szekeres metrics with cosmological constant 

Since H^^ (^) ^ - 0{6p, S), we have that: R-'^nens 0{Sp+'' , 3^+") whereas R-^kEs 0{Sp, S). Thus, the non- 
spherical perturbation to the energy density is always of sub-leading order compared to the first order in spherical 
perturbation. The first-order, non-spherical, metric perturbation appears at 0{S'^); however, this is equivalent to 
a coordinate transform on (r, 9, ip) and does not source a non-spherically symmetric physical perturbation to the 
dilaton evolution at this order. For the dilaton field, </>, then, the first non-spherical perturbation is always sourced 
at subleading order compared to the first spherically symmetric one. 

In the exterior we can apply the boundary condition that cj) ^ (j)c{t) as r oo. In addition to this we also have the 
stronger condition that, as S —^ 0, the inhomogencity should disappear and (j) — > 4>c{t)- Thus to zeroth order in the 
exterior (f> ~ 0c (0 + 0{dP, 6). Since we are only really interested in the behaviour of (/> in the interior we do not need 
to calculate the higher-order terms in the exterior limit explicitly, we only need to know enough about their behaviour 
to be able to perform the matching in some intermediate scaling region. We will consider that behaviour later. 

B. Interior Limit 

To lowest order in the interior region, we write Z ~ S'^RsfJ.i, and :— 5'' Rj^bi, where i = {x,y,z}; q and q' are 
determined by specify a particular matter distribution. The condition that k£ > everywhere requires q' > q. From 
the exact solutions we find: 

k{r) ^ j2/3^^ (1 ^ ^ „ (^9)) ^ 1^^5/3^ ^ 

where ko{6T) = {2m/Rs) {-n/iHoh -h 5T)f'^. In refs.[lll we sought to remove the effect of the 0{5'^l^) in k{r) by 
a transformation of the time coordinate. However, it is not clear that this new time coordinate is well-defined near 
a black-hole horizon; we now believe this procedure to have been technically incorrect (although it did not effect 
the results). We correct it in this work by implementing the 0{5^/^) correction differently, li q' > q then to the 
next-to-leading order, we need only consider the spherically-symmetric modes to find interior expansion of 0. We 
could also include a non-spherical vacuum component for (j) at next-to-leading order; however, this will be entirely 
determined by a boundary condition on R = Rs and the need that it should vanish for large R. To find the leading- 
order behaviour of the (j)^T we need to know ip at next-to-leading order. Hence, the only case where we must explicitly 
consider non-spherically symmetric effects is when q' — q, i.e. KSns = O^kEs)- In what follows it is natural to consider 
the spherically symmetric and non-spherically symmetric modes of (j) separately. 

Before we can solve the equations in the interior limit, we need a boundary condition at R — Rg. At leading 
order we take this to be: 

Rl (l - dR<j>o\n^R^ - 2mF {^o) 
dR'R'^B,^{(l)Q{R', t))Ke{R'), 

The no-hair theorem for black holes implies that F (0o) = 0, however for bodies where 2m/ Rg ^ 1 we expect 
F (0o) ~ B^^{(j)c)- At higher orders we find the flux F by perturbing the above expression as explained in 11]. The 
zeroth-order mode of (f) in the interior is then found to be: 

The matching procedure gives 4>e (ST) — (j)c{t)- 

1. Spherically Symmetric Perturbations 

In we considered the spherically symmetric perturbations of (p that occur at order 5* and at order 5. Here, 
we will also consider the perturbations at order 6"^^^. The perturbation at order 5 is sourced by the □o0e(i5T') term. 
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where Dq is the d'Alembertian of the Schwarzschild metric. This effect of perturbation acts like a drag term, and is 
equivalent to the change (j)e{ST) (p^^ST') with: 



T' = T - 2 



2m 

i?7 



3/2 



2m 



-In 



2to 



where C is some constant which is determined by the matching process. When the matching is performed one finds: 



T' r^T- 



^2(1 - ki^'^t) - {2m/Rs)y^i\{i',T)) 



where £,,t{£,,T) — R^t{R = Rs(.,t = io + LjT). When the 'star' is actually a black hole, the above expression 
reproduces, to leading order in 5, the result found by Jacobson j^^. As for the 0{5'^/^) correction to the metric 
coming from fc(r), we find that it only effects at order 0{6^/^ ^ (5''+^/^); as such, they can be ignored since they are 
always smaller than the effects that we have included. 

In we calculated the order 5'^ correction to t/)'-*^ (T, ^) . In cases where the surface of our 'star' is far outside 

its Schwarzschild horizon (2m/i?s), (j)f^ is given by: 



2m 

lis 



B( 



DiT) 




In the cases where 2m/ R^ 1 it was not possible to find a closed analytical expression for 0^'''; however, from the 

its equation of motion, it can be easily seen that will be of the same order of magnitude as the above expression. 
The function D{T) is a constant of integration, and it must be found via the matching procedure. Before we perform 
this matching, we will consider the non-spherically symmetric modes. 



2. N on- Spherically Symmetric Perturbations 

Non-spherically symmetric modes in the interior approximation to will be sourced at order 5'^ , where q' > q. 
We studied these modes and found that, to order (5*, they only have a dipole moment. We write the order (5* 
non-spherically symmetric, modes as (5* where: 

(f)\^"^ ■= 0^'')^(^, T) cosd + (t>\^'^''{£„ T) sin 6* cos if + (t>^p^{£„ T) sin 6* sin if. 
As before, when 2m/ Rg ^ 1, we can find analytic expressions for these modes: 

—B^^ (0?) e j dr^'^ + — 5 ^ (0?) J, 

- ^Pi^o)^ r dv'k{v')e + ^ + D,{T)^ + 0{{2m/R,n 

where i — x, y z. When 2m/ Rg « 1 the above expression can be seen as an order of magnitude estimate for the 
^{q rpj^g Di{T) and Ci{T) are constants of integration, with Di[T) determined by the matching procedure. The 
value of Ci{T) should be set by a boundary condition on R = Rg. We cannot specify Ci{T) exactly without further 
information about the interior of our 'star' in R < Rg. If we assume that the prescription for the sub- leading order 
boundary condition given above is correct then we find: 



i^')' ^ (0?) e / d^-^ + , (0?) ± drj'eb^w) (8) 



o Aq')i\ 2m h 



F{4>o)+O{{2m,/Rgf 



Rs r;i/2 



From now on, we set Ci = 0, for simplicity. Even when this is not exactly satisfied, we do not expect the magnitude 
of Ci or d^T to be larger than any of the other terms in 0^^^' or respectively. 
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C. Validity of Matching Procedure 

Before we can apply the matching procedure, we must ensure that it is appUcable to our problem i.e. that there 
exists some intermediate region where both the interior and exterior approximations are simultaneously valid. We 
considered these conditions in ref.^. We define coefficients n > and > by /ii (ry) ^ 77" and bi (ry) cx rj'^^ as 
7] ^ 00 respectively. Writing H^^Xi ~ 5P'li{p), we also define coefficients m and fi by zi{p) ^ p™ and li{p) cx p^^' 
as p — > 0. For both the exterior and interior to be simultaneously valid in some intermediate scaling region where 
ViXii ^ (5^" with < a < 1, we need there to exist some a such that: 

fn^ P \ ^ ^ 1 

max U, 1 < a < — 

\ 1 — 771 y n 

max(p^ + (1 — ce){.fi + ™)) > — J3 if P < 1, 

i 

max(p,- + (f - a) ft) > -f if p > 1. 

i 

a — max q/di > 

i 

These conditions can, in almost all cases, be rephrased as: 

limi^^KAe ^ off) (9) 

(5^0 

lim2(TO + Z)/i? = 0(1) (fO) 
as 5 ^ 0, with LfL]f"{t - to),LfL]^-°'R held fixed, for all a G (0, f ). 



D. Matching and Results 

We are interested in time derivatives of the </> field. From the expression for (p in the interior, we see that t' = LjT' 
seems to play the role of a natural time coordinate. At radii where 2m/ R <C f , the interior metric is close to diagonal 
when written in {f , R) coordinates; it is in this sense a natural time coordinate for an observer at fixed R. In this 
region, t' coincides with the standard Schwarzschild time coordinate. As R ^ 00, t' t. We therefore consider ip^t'- 
Whenever the required conditions of the previous section hold, the matching procedure is valid, and wc find 

4,,t' {r, t) « 4"], + <5VS' + ^'^VK'^ + o{5\ 

in the interior, where the « sign means that whilst this is not a formal asymptotic series (since there may be excluded 
terms that are bigger than some of the included ones) this is a good numerical estimate since at least one of the 
included terms will he bigger than all the excluded ones. The terms in this expression are given by: 

4°t' ~ ^eW + At(r,O0c(i), (11) 



. , . / w MRt)R^ - RlMRt)\R=R , , 



where Ai?,f = - BR, 



8H%--BA^.)\ \ dr'ARAR,.e)ir',t) + ( 1 - (^-^-^^^^^^-^^ ) + ^ ^''^ 




and 



SU\''J - -^5,0 (0?) R rdr'R,rR,t^^^^^p^ - l:B,^ (0?) dr'i?,.i?,,i?2«£„,(r', t) 

+ i^o) ^2 / dr'R.rR^tR'^KSnsir', t) - (0o) RR.tne^sir, t) + O 
•J Jr-r„ o 



(14) 




Equations (|f 3|l and (|f4|l are derived for the case 2m/ Rs ^ 1. For 2m/ Rs ~ 1, they are accurate when R ^ 2r7i and 
otherwise provide order-of-magnitude estimates for S'^(f)''flt and S"^ (ji'it' respectively. 
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We can now evaluate these terms to find out when (j){x,t) « (j)c{t) and also state the precision to which this 
approximate equality holds. In many cases, however, a lot of the terms in the above expression are negligible or 
cancel, and so we can find a more succinct necessary and sufficient condition for 4){x,t) w ipcit) to hold. When 
2m/ Rs > 1 we expect F (0o) ~ ^,0 {4>fj (1 + 0{2m/Rs)), and so: 

(l>i,t' -cl)c,t ~ ~B^^(c^,) [ dr'R,rRtKAes{r',t)~lB,^{^,)R f dr' R^rR.t^^^^^^^^ (15) 

- ^B,4,{(j)c)RR,tKens{r,t) + At{r,t)4,,{t). (16) 

We will refer to this last term as the drag term, and it is responsible for the local value of 4> lagging slig htly behind 
the cosmological one - this effect was first observed by Jacobson in the study of gravitational memory, [2J| . Whenever 
the cosmological (p is not potential dominated, and our 'star' resides in a local overdensity of matter, the drag term 
will be negligible compared to the other terms in this expression. If the potential term dominates the cosmological 
evolution then it is possible for the drag term to give the dominant effect, even if we have an local over-density of 
matter. But whenever this happens we always have |A<</)c/0c| 1 and so we have (j){x,t) « 4>c{i)- 

Independent of the nature of overdensity, we saw in our previous papers that potential domination of the cosmo- 
logical 4> evolution acts only to strengthen the degree to which (j){x,t) « (pdt). 

Even if we ignore the drag term, the above expression for (jjj^t' — 4'c,t is still rather unwieldy. In almost all cases, 
integrating over the non-spherically symmetric modes of ke in the same way as we did for the spherically symmetric 
ones only acts to increase l^/,*' — </>c,t|- Therefore, we define the quantity J by: 

I := B,(0.) /^di^. "^ax,,,(singA(.s)) ^ ^^^^ 



where v is the radial velocity of the matter particles (i.e. v = i?.t). When the background spacetime is Szekeres-Szafron 
we have that: 



with equality in the spherically symmetric case. The strong inequality X ^ 1 is therefore a sufficient condition for 
(j){x, t) « ipcit), and the value of I gives a measure of the amount by which (p^t' and (pct differ. In addition to having 
shown this for cases where the background is Szekeres-Szafron, and that matching conditions hold, we also conjecture 
that, even if the matching conditions formally fail, and for other classes of spacetime, that I <C 1 is a sufficient 
condition for 4i{x,t) « (pcit)- 

VI. COLLAPSING CASE 

When a spacetime undergoes gravitational collapse it is possible for a black-hole to form inside the collapsing region. 
In the Tolman-Bondi model a black-hole horizon appears when 2M (r) /R = 1. If we have kAeR^ = {ne — k£c)R^ <C 1 
outside the horizon then we can apply the results of the previous section, taking the surface of our 'star' to be the 
black-hole horizon. This is also true for any virialised region in the interior of the collapsing region, not just for black 
holes. 

The results of the previous section can also be extended to the case where the collapsing interior region has no 
central black-hole or virialised region. For simplicity we consider only spherically-symmetric, dust-plus- A cosmologies 
i.e. Tolman-Bondi models. We do not require the big bang to be simultaneous for all observers. This extension 
requires that curvature of the interior spacetime be in some sense weak so that the metric is Minkowski to zcroth 
order. We require: 

R^KAe<.l, 2AM/R<.1 

everywhere; Ae = e — Cc and AM = M— ^kEc- in the interior. In the exterior we assume, as before, that the spacetime 
is FRW to zeroth order. It is clear that in this model the following parameters, Si and ^2 are everywhere small: 

2M 

di{R,t) = i?2(Ke-KeJ = — ^-3^^„i^2i?^ 

R r 

6,{R,t) = ^-n^H'R\ 
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In addition, the following parameter, 5^, is small in interior region but 0(1) in the exterior: 

In the interior 63 <C Si, 52- For the purposes of our asymptotic expansions we treat and and 62 as being of the same 
order. For the interior to be collapsing wc need k{r) > 0. The condition that R\ = —k{r) + ^2 + {^m + f^A)^3 > 
implies that k{r) < 0(^2), and R\ ~ 0(62) in the interior. 

We will perform the matching in an intermediate region where 61 S2 S3 <^ 1. It is clear that with these 
definitions that such an intermediate region must always exist. 



A. The Interior 

In the interior we write the metric as: 

,2 {l-S2+Ss{nm + ^x))dt^ , 2R,tdRdt dR^ r2 r ,.2 ^ • 2 . , 2. 

which is flat spacetime to lowest order in the Si. The dilaton, cj), obeys: 

-R^n^ = B^4cj))iSi + SQmSs) + V4(l>)R^, 

and, in line with our previous assumptions, we have V^^{(t))R^ ~ (^(nA^a)- We note that dtSi ~ 0{S\^'^S2) = o{S\^^) 
and so Si is quasi-static; as such we expect </> to also be quasi-static in the interior. We can solve the equations for cj) 
order-by-order in the interior, requiring (as a boundary condition) that (j) is regular at R = Q: 

rR ,1 D/ rR 



S) + B,^{^e) '^S2{R',t)+Mt) dR'{R,t-HR') (18) 
+ ^ [i?V^(</)e) + 3B ^0^53 + {k{t) + 3F0e(i))i?'] + O (si, Si, Si, {R^,f, S^iR^S) , 



where 0e(i) is C(l) but quasi-static i.e. i?(/>e = o{S-{ , 82 ); the third term is 0{52 R4>e)- Since the above expression 
is not a formal asymptotic expansion as such we cannot be sure that the neglected terms are smaller than all of the 
included terms; indeed wc shall see that the matching ensures the vanishing of the term in [..]. This is because we do 
not know precisely how the sizes of 83 and R4>e relate to those of 5i and 82- What we do know is that, in the interior, 
the excluded terms are smaller than at least one of the included terms. The limits C and D as well as 4'e{t) must be 
found matching the interior expansion to the exterior one. 



B. The Exterior 



In the exterior we define a coordinate q = R/a{t) r where a{t) is the FRW scalar factor. In {t, q) coordinates the 
metric reads: 

ds^ = dr (1 + 0{Si, {AkY) + — — - - a^£)^{dr + sm edcj)^} 

1 — k(r) 1 — k[r) 

where A/c = k{r)-koe^, and ko = \imr^^ k{r)/r'^, (l-n„-0A)i?2 = -k^/a^ and 2 - HR) ~ {-Ak + S2)/HR. 
As R,g ^ 00 and the inhomogcncity is removed (i.e. Ak, 61,82 0) we require that 4> (pcit)- As with the 
virialised case our only interest in the sublcading order behaviour of (j) in the exterior is so as to match it to the 
interior approximation. It is only necessary therefore to consider how the exterior approximation to behaves in the 
intermediate region where all the Si are small. In the intermediate region the exterior approximation is: 

<!> ~ .^e(i) + (I>s{q, t) + B ^(./-c) '^S2{R', t) + 4>c{t) r dR'{R,t - HR') + 0{8lSl, Sl), 
where (/>s(£», t) — o(l) is some vacuum mode i.e. \^FRW<t^s = 0- 
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C. Matching and Results 

Now we match the interior and exterior expansions and find that C = D = oo, (j)e{t) — (f>c(t) and (ps = 0. The 
requirement that (t>eit) = (t>c{t), combined with the cosmological evolution equation for (jjcit) ensure that the term in 
[..] in eqn. (|18|l vanishes. We have found that the matching interior approximation is therefore given by: 

~ Mt' = < + At) + S ^S^iR', t) + 0{5l5l 51). 

where the lag in the time coordinate, At, is given by: 

At^ AR'{R,t- HR'). 

J oc 

This coincides (to leading order) with the expression for the virialised case, eqn. (|12|l . for Rs = 0. It seems natural, 
in the interior, to consider the time derivative of w.r.t. t' . As noted in the previous section, t' will look like the 
Schwarzschild time coordinate near the surface of a massive body (that is far outside its out Schwarzschild radius), 
and t' ^ t as i? — > oo. We find that: 



B. 

(j),t'{r, t') - (bc,t 0(0c) / di?'A(i?,tKe)(i?', t) + 4>,{t)At + 0(61 6lSl) 

J oc 

where A(i?_tK£)(i?, t) = R^ti^siRit) — HRkEc- We note that this is the same as the spherically-symmetric limit of the 
result found in eqn. (|15|l for the virialised case. 

This completes the extension of our analysis to the case of spherically symmetric collapsing spacetimes. It is clear 
that the quantity I, defined in the analysis of the virialised case, will be also be a good measure of |(0 — (/'c)/</'c| in 
the collapsing case. The validity of the matching procedure is this case is assured by the condition that (5i,(52 ^ 1 
holds everywhere. 

VII. RESULTS AND CONSEQUENCES 

We now consider the astronomical consequences of our results for observations here on Earth, and answer the basic 
question of whether local experiments will detect cosmologically varying constants. We can evaluate the quantity T 
explicitly for an Earth-based experiment assuming the varying constant to be the Newtonian gravitation "constant" 
G governed by Brans-Dicke theory (since in this case the cosmological evolution of is easy to solve). We expect 
similar values for BSBM, BM and other non-potential-dominated theories for varying a and /i 

We will consider a star (and associated planetary system) inside a galaxy that is itself embedded in a large galactic 
cluster. The cluster is assumed to have virialised and be of size Rdust- Close to the edge of the cluster we allow for 
some dust to be unvirialised and still undergoing collapse. There are three main contributions to I coming from the 
star, the galaxy, and the galaxy cluster, respectively, and of these the galaxy cluster contribution is by far the biggest, 
this can be understood by noting that the galaxy cluster is the deepest gravitational potential well, and the galaxy 
and star are only small perturbations to it. The contribution to T from the galaxy cluster is found to be: 



^ciust ^ ^Ho{s — l/2) ^\/2MciustRc 



~^clust 



iust ' 



_ m/{'2s~l)]vl^A^ + z.„f'^Al{l _ 1 ni X in-5 1^ [3/(2.-l)](l + z..,,)5/2Al^ \ ( hM^V'' 

= 4.95 X 10-4[3/(2. - l)]f77„V2 (i + ,_)3/2^ 

« 1.61 X 10-3[3/(2s - l)]f7„i/2(i ^ z„^)3/2 ^ 1 

where we have used ?>Mciust/^Rciust = vliust = ^CTu and nedust = &Mciust/ Rliusu is the 1-D velocity dispersion 
and Avir ~ 178 is the density contrast between the cluster and the background at virialisation. In the final line of the 
approximation we have used the representative value ay = 1040kms~^ Vyir = IBOOkms"^ appropriate for a rich 
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cluster like Coma, [26j- Taking a cosmological density parameter equal to ilm = 0.27, in accordance with WMAP, we 
expect that for a typical cluster which virialised at a redshift Zy^r 1, we would have Tciust ~ 0.31[3/ (2s — 1)] x 10~^. 
The term in [..] is unity when s — 2, i.e. 2GM/R const; such a matter distribution is characteristic of dark matter 
halos. Different choices of s > 1/2 only change this estimate by a factor that is 0(1). We note that, since 2GAM/R 
(where AM — M — ^ScR^) is required to be small as i? — *■ oo by the matching conditions, the model used here is 
only valid for s < 2 and hence the singularity in Tdust at s = 1/2 is fictitious. If we were to have Tdust ^ 1 then we 
would require a large virial velocity: Vdust ^ 32, 400[3/(2s — 1)]~^/^(1 + z^ir)~^^'*kms~^. 

It is clear that in theories like Brans-Dicke, which have their cosmological evolution dominated by the matter-to- 
dilaton coupling, _B 0, the local time variation of </) and the associated constant differs from its cosmological value by 
at most about 1%. In theories where the potential dominates the cosmological evolution this result becomes even 
stronger and we expect any deviations to occur only at the QA\B^^{(f)c)/V^fi,{4>c)\% level, where |B_0((/)c)/V^0(0c)| <C 1. 

We have also seen that I is also good measure of \{(t) — (\)^l4>c\ inside spherically-symmetric regions that are still 
undergoing gravitational collapse. Let us now evaluate T for the case of a collapsing cluster in Brans-Dicke theory in 
the matter era. We assume that the cluster is approximately homogeneous. We further assume that when the cluster 
eventually virialises, at time tmr-, it has a virialisation velocity Vmr- We use the spherical approximation detailed in 
Chapter 5 of to model the collapse of the the cluster. We define by t = tyir{6 — sin0)/27r, and find: 

where A^^^ = IStt^ « 178 is the density contrast at virialisation when 9 ~ 2tt. When 9 < 3Tr/2 f{9) — sin6'(l — cos 9)~^ 
and g{9) = sin6'; for 9 > 37r/2, f{9) = g{9) = —1. In this evaluation we have included the effect of the 'drag term'; 
this is important up to turnaround but it becomes negligible soon afterwards. Turnaround occurs aX 9 — T:,t — ttum- 



lit) 



5vt 



.1/2 



f {9)(9- Sin 9) ^ V2ng(9) 
'2V27T AV^(6'-sin 




FIG. 1: Plot of {cp^t' ~ 4>c,t)/4>c,t vs. time for Brans-Dicke theory at the centre of a collapsing cluster with = 1800km s ^ . 

In deriving the above expression we have used SMdust/^Rvir = 'L'vir^ where Rvir is the radius of the cluster after 
virialisation and Mdust its mass. The conditions required for the matching procedure to be valid are equivalent to 
10w^j^/(l — cos 6*) ^ 1, and it is clear that this will not be satisfied all the way down to 6* = 0. For Vmr = 1800kms~i, 
our method will be valid for 9 > 0.027 and for the matching conditions to hold from turnaround to virialisation 
we require Vmr ^ 95000km s^^. Assuming that the cluster virialises at an epoch that is close to the present day, 
this bound on Vyir translates to requiring Ryir <C 432/i~i(l -f- Zyir)~^^^M-pc, where Hq = 100/ikms~i Mpc~^. We 
observe that X is small up until turnaround and then grows quickly until virialisation. At turnaround X = 0, and at 
virialisation we find 
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For the final evaluation we have taken Vmr = ISOOkms ^ (as appropriate for the Coma cluster). The vanishing of 
X at turnaround is specific to Brans-Dicke theory, more generally: I{tturn) = 40?;^j^[^c — -B,0((/'c)«ec/-ff'/>c]/277r^. In 
theories where the matter coupling is strongly dominant cosmologically, \B'^^{4>c)Ktc\ ^ IK<>('/'c)|) we find I{tturn) ~ 
16O«,2,,|B^^(0e)|/277r2 « 1. ' ' 

Our results differ greatly from those that were found using the spherical collapse model used in [23; where I{tmr) ~ 
200. In that model, the spatial derivatives of (j) were assumed to be negligible and are neglected. However, this is 
can only be a realistic approximation when the collapsing region is as large as the cosmological horizon; for a cluster 
virialising today that would require R^r ^ 5Gpc. Since our method will fail for R^r ^ 432/i~^(l + z^ir)~^^^Mpc, 
there must be some region of intermediate behavior, SOOMpc < R^r S 5Gpc, that is not described by either the 
spherical collapse model or our present analysis. We derived these results for Brans-Dicke theory, where (p oc G~^, 
however we should expect similar numbers for all varying-constant theories where the cosmological dilaton evolution 
is dominated by its matter coupling, B ^pHEc- In potential-dominated theories, the above numbers will be reduced 
by a factor of \B,ci,{(f>c)K£c/Vcii{(f>c)\ ^ 1- As in the post-virialisation case, potential domination of the cosmological 
evolution strengthens the amount to which local experiments will see cosmologically varying constants. 

In conclusion: we have used the method of matched asymptotic expansions to find a sufficient condition for the 
time-variation of a scalar field, and any related varying physical 'constants' whose variation is driven by such a field, 
to track its cosmological evolution. We have extended our previous analyses by allowing 'local' also to include being 
inside some spherically-symmetric collapsing region. We have also proposed a generalisation of our earlier condition 
for local variations to follow global cosmological variations. We conjecture that this new condition is applicable to 
scenarios more general than those we have explicitly considered here. We have seen that this sufficient condition 
is always satisfied for typical distributions of matter, and we have provided a proof of what was previously merely 
assumed: terrestrial and solar system based observations can legitimately be used to constrain the cosmological time 
variation of many supposed 'constants' of Nature. 
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